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ABSTRACT 

We present high spatial resolution optical and near-infrared imaging obtained using the ACS, 
WFPC2 and NICMOS cameras aboard the Hubble Space Telescope of 31 24/im-bright z « 2 Dust 
Obscured Galaxies (DOGs) identified in the Bootes Field of the NOAO Deep Wide-Field Survey. 
Although this subset of DOGs have mid-IR spectral energy distributions dominated by a power-law 
component suggestive of an AGN, all but one of the galaxies are spatially extended and not domi- 
nated by an unresolved component at rest-frame UV or optical wavelengths. The observed V — H 
and I — H colors of the extended components are 0.2 — 3 magnitudes redder than normal star-forming 
galaxies. All but 1 have axial ratios > 0.3, making it unlikely that DOGs are composed of an edge-on 
star-forming disk. We model the spatially extended component of the surface brightness distributions 
of the DOGs with a Sersic profile and find effective radii of 1 — 6 kpc. This sample of DOGs is 
smaller than most sub-millimeter galaxies (SMGs), but larger than quiescent high-redshift galaxies. 
Non-parametric measures (Gini and M20) of DOG morphologies suggest that these galaxies are more 
dynamically relaxed than local ULIRGs. We estimate lower limits to the stellar masses of DOGs 
based on the rest-frame optical photometry and find that these range from ~ 10^~^^ Mq. If major 
mergers are the progenitors of DOGs, then these observations suggest that DOGs may represent a 
post-merger evolutionary stage. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: high-redshift 



1. INTRODUCTION 

One of the most important questions concerning the 
evolution of galaxies is when and how the most massive 
galaxies formed. It has been known since the analysis 
of the InfraRed Astronomical Satellite (IRAS) data that 
in the local universe the most bolometrically luminous 
galaxies have their spectral energy distributions (S EDs) 
dominated by infrared (IR) light fSo ifer et al.lll986[ ). sug- 
gesting that these systems are highly obscured by dust, 
which absorbs ultra-violet (UV) and optical light and 
re-radiates it in the IR. While these ultra-luminous IR 
galaxies (ULIRGs) are rare in the local universe, they be- 
come an i ncreasingly important phenomenon at high red- 
shift (e.g..lFranceschini et al.l l2001l:l Le Floc'h et ai.l l2005l: 
iPerez-Gonzalez et al.ll2005f ). 

Following the launch of the Spitzer Space Tele- 
scope, numerous investigators have identified and stud- 
ied populations of high redshift galaxies that are IR- 
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bright yet optically faint (" Yan et ani2004l: iHouc k et al 



2005|:IWeedman et a l. 2006b : iFiore et al.l l2008: De v et al 



20081) In particular. Dev et all (|2008l ) and iFiore et al 
()2008l ) present a simple and economical method for 
selecting these systems using only i?-band and 24^m 
Multiband Imaging Photor neter for S pitzer (MIPS; 
iRieke et al.l l2004f ) data. iDev et all (j2008l ) employ 
a color cut of i? — [24] > 14 (Vega magnitudes; 
wF,,(24^m)/i^i.(i?) > 1000) to identify objects they caU 
Dust Obscured Galaxies (DOGs) in the Bootes field of 
the NOAO Deep Wide-Field Survey (NDWFS^^; Jannuzi 
et al., in prep.; Dey et al., in prep.). 

The broadband photometry, redshift distribution, and 
number density of the DOGs imply that they are under- 
going an extremely luminous, short-lived phase of stellar 
bulge and nuclear black hole growth and may be the pro- 
genitors of the most luminous (~4L*) present-day galax- 
ies. Ground-based photometry from the NDWFS sug- 
gests magnitudes of i? w 24 - 27, / w 24 - 26, and K k. 
17.5 - 20.5 for the sample of DOGs with F^{2Afim) > 
0.3 mJy. DOGs are relatively rare, with a surface den- 
sity of «0.089 arcmin^^ for sources with i^24 > 0.3 mJy. 
Spectroscopic redshifts determined for a sub-sample of 
DOGs using the Deep Im aging Multi-Object Spectro- 
graph CDEIMOS: [Faber et al.li2003[) a nd the Low Res- 
olution Imaging Spectrometer (LRIS; lOke et all Il995f) 
on the telescopes of the W. M. Keck Observatory (43 
DOGs), as well as the Infrared Spectrometer (IRS; 
IHouck et ai] I2004D on Spitzer (43 DOGs) have shown 
that the DOGs have a redshift distribution centered on 
z « 2 with a dispersion of <Jz ~ 0.5. While DOGs are 
rare, they are sufficiently luminous that they contribute 
up to one-quarter of the total IR luminosity density from 

^'^ http://www.noao.edu/noaodeep 
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all z ^ 2 galaxies, and constitute the bulk of ULIRGs at 
z ~ 2 iDe^etaDllQQl). 

Based on their observed properties, iDev et all ()2008D 
suggest DOGs may represent a transition stage be- 
tween sub-millimeter-selected galaxies (SMGs) and un- 
obscured quasars or galaxies. Evidence in support of 
this scenario is that DOGs and SM Gs have similar spac e 
densities and clustering properties (jBrodwin et al.|[2008l ). 
An important test of this scenario is to study their mor- 
phologies with high spatial resolution imaging. For ex- 
ample, one of the primary motivations for the merger- 
driven scenario for the formation of ULIR Gs is their dis- 
turbed structure at optical wavelengths (jSanders et al.l 
[T988a) . Studies of numerical simulations of galaxy 
mergers have suggested that they can produce very 
red, luminous syste ms that a re highly dust-obs cured 
(jJonsson et al.ll2006l) . Recentlv. iLotz et al.l (|2008b[ ) have 
applied non-parametric methods of quantifying galaxy 
morphologies to similar merger simulations and have 
found that mergers are most easily identified during the 
first pass and at the final coalescence of their nuclei. 

In addition to identifying merger activity, morpholog- 
ical information can constrain the size-scale of the emit- 
ting region. Sources with active star formation on sev- 
eral kiloparsec (kpc) scales have larger sizes than ob- 
jects dominated by an Active Galactic Nucleus (AGN) 
or a very compact, nuclear starburst. Studies of Dis- 
tant Red Galaxies (DRGs) have shown a relation between 
star formation and size at rest-frame optical wavelengths, 
in the sense that quiescent DRGs are all very compact 
with effective radii (i?off) less than 1 kpc, while active 
DRGs tend to be more extended (1 < i?off < 10 kpc; 
iZirm et al. I [20071 : 1 Toft et al.l[2?)nl . Analysis of SMGs in 
GOODS-N shows extended emission on scales of 5-15 kpc 
(|Pope et al.ll2005f) . Rece nt NICMOS i maging of a sample 
of 33 high-z ULIRGs bv lDasvra et al . (2008) has shown 
these extreme objects (which are similar in their selec- 
tion criteria to DOGs) to have effective radii in the range 
~1.5-5 kpc. About half of their sample shows signs of 
interactions, but only 2 are merging binaries with a lu- 
minosity ratio <3:1, i.e., qualifying as major mergers. 

High spatial resolution imaging of the DOGs is es- 
sential to understanding their relation to other galaxy 
populations as well as their role in galaxy evolution in 
general. We have begun an effort to obtain high reso- 
lution imaging using laser guide star and natural guide 
star adaptive optics on the Keck telescopes. These (on- 
going) efforts have resulted in high resolution iiT-band 
images of a handful of DOGs found near bright stars 
(|Melbourne et all l2008l , in prep.). A complementary 
method of obtaining deep, high spatial resolution imag- 
ing is with the Hubble Space Telescope (HST). With the 
Advanced Camera for Surveys (ACS) and the Wide Field 
Planetary Camera 2 (WFPC2), we can probe the rest- 
frame UV emission of the DOGs that is sensitive to the 
ionizing sources associated with on-going star formation. 
Meanwhile, NICMOS data allow the study of the rest- 
frame optical morphology, which better traces the stellar 
mass and dust-enshrouded AGN. 

In this paper we present ACS/WFPC2 and NICMOS 
images of 31 DOGs and analyze their morphologies. 
The DOGs studied in this paper have spectroscopic red- 
shifts from either Spitzer/IRS DEIMOS/LRIS, were se- 
lected primarily based on their large 24/im flux densities 



(F24/im > O.SmJy), and have power-law SEDs in the mid- 
IR. In a future paper, we will study a sample of DOGs 
with fainter 24/im flux densities that have mid-IR bump 
SEDs (Bussmann et al., in prep.). In section[2]we detail 
the sample selection, observations, and data reduction. 
Section [3] contains a description of the methods we use 
in our morphological analysis, and in section 31 we re- 
port the results this analysis. In section [SJ we estimate 
some intrinsic properties of the DOGs in our sample and 
we compare our findings with what is seen in other high 
redshift galaxy populations. Finally, we present our con- 
clusions in section [6l 

Throughout this paper we assume 
Ho =70 km s-i Mpc'^ = 0.3, and fix = 0.7. 
At z = 2, this results in 8.37 kpc/". 

2. DATA 

In this section, we describe our sample selection and 
give details regarding the i/ST observations and our data 
reduction procedure, as well how we measure our pho- 
tometry. Finally, we show postage stamp images and 
provide a brief qualitative description of each target. 

2.1. Sample Selection 

As outlined in sect ion [U a sample of Ri2600 DOGs 
from lDev et al.l ()2008[ ) was originally identified using the 
9.3 deg2 Bootes Field of the NDWFS. For details of 
the selection crit eria and photome tric analysis, we re- 
fer the reader to iDev et al.l (|2008D . In this paper, we 
analyze HST imaging from program IIST-GO10890 of 
31 of the brightest DOGs at 24/iTO (all have F24pm > 
0.8 mJy). The bolometric luminosity of DOGs with 
bright 24/im flux densities is typically dominated by 
AGN emission, while the opposite is true for 24/tm faint 
DOGs (0.1 mJy < i^24u m < 0-3 m Jy), which arc domi- 
nated by star-formation ijPope et a l. 2008). Additionally, 
IRAC photometry shows that the objects in this paper 
are dominated by a power-law component in the mid-IR. 
The most likely cause of this er riission is the presen ce of 
warm dust heated by an AGN (|Donlev et al.ll2007[ ). 

Shallow X-ray coverage of the Bootes field ex- 
ists and has yielded a full catalo g of X -ray sources 
(iMurrav et al.l 120051: iKenter et all 120051: iBrand et al] 
120061). Within a 2" search radius, two of the DOGs 
studied in this paper (SST24 J143102. 2-^325152 and 
SST24 J143644.2-f 350627) have a single X-ray coun- 
terpart, and one DOG has two counterparts (SST24 
J142644.3-h333051). A full analysis of the X-ray data 
is beyond the scope of this paper, but these basic re- 
sults suggest that most DOGs are either not strong X- 
ray emitters or are heavily obscured. The latter view is 
supported both by mid-IR spectral features and the fact 
that this subset of 24/im bright DOGs shows some of the 
reddest R — [24] colors of the entire DOG population. 
Figure [1] shows the color-magnitude diagram in i? — [24] 
vs. [24] space for the full DOG population in Bootes and 
highlights the subsample of objects studied in this paper. 

Previous work has shown that objects dominated by 
a power-law signature in the mid-IR tend to have AGN 
indicators in their mid-IR spectra, usually silicate ab- 
sorption but no PAH emission (W ccdman et al. 2006a!; 
iPolletta et all 120081 : iBrand et all l2008f ). Indeed, IRS 
spectra of these sources have revealed redshifts based 
on the 9.7/tm Silicate absorption feature, and all are 
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Fig. 1. — R — [24] color vs. 24/xm magnitude distribution for 
DOGs in the NDWFS Bootes field. Bottom and top abscissae show 
the 24/xm magnitude and flux density, respectively, and the left and 
right ordinates show the color in magnitudes and the F2Af_im/ Fji 
flux density ratio, respectively. Black dots and upward arrows show 
the full sample of DOGs, with and without an i?,-band detection, 
respectively. The subsample studied in this paper is represented 
by red circles (open symbols show lower limits). 




Fig. 2. — Distribution of redshifts for DOGs in the Bootes Field 
with spectroscopic redshifts (either from Spitzer/IRS or Keck 
DEIMOS/LRIS). The redshift distribution of the sub-sample of 
objects studied in this paper is shown with the hatched histograms 
and is representative of the full sample of DOGs. 

located at z ~ 2. Of t he 31 object s in t his sam- 
ple, 17 have spectra from [H ouck ct al.l ()2005D . 2 have 
spectra from IWeedman etal- (,2006bh . and the remain- 
ing spectra will be presented in future wor k (Higdon et 
al. in prep. ). Subsequent Keck/NIRSPEC (|Brand et all 
I2007f ). Keck/LRIS, and Keck/DEIMOS spectroscopy has 
yielded more precise redshifts for 4 of the DOGs. The 
redshift distribution of the sample studied in this pa- 
per compared to the overall distribution of spectroscopic 
redshifts for the DOGs from the Bootes field is shown in 
Figure [H 

2.2. Observations 

The 31 DOGs we study here were observed with HST 
from 2006 November to 2008 February. Nine were im- 
aged in the Wide Field Channel (WFC) mode of ACS 
(|Ford et al.l Il998l ) before the failure of the instrument. 



We have observed th e remaining 22 DOGs with WFPC2 
(|Trauger et al.lll994D . All 31 DOGs were observed with 
the NICMOS NIC2 camera. Table [U summarizes the de- 
tails of the observations. All data were processed using 
IRAF^^. In the following sections we provide more details 
about the processing of the ACS, WFPC2, and NICMOS 
images used in this paper. 

2.2.1. ACS 

Each DOG was observed over a single orbit through 
the F814W filter using a four point dither pattern (ACS- 
WFC-DITHER-BOX) with a point spacing of 0.265", a 
line spacing of 0.187" and a pattern orientation of 20.67°. 
Total exposure time was «2000 sec. Bias-subtraction 
and flat-fielding was performed using the standard ACS 
pipeline (Pavlovsky et al. 2004). The MultiDrizzle rou- 
tine was used to correct for geometric distortions, per- 
form sky-subtraction, image registration, cos mic ray re- 
jectio n and final drizzle combination (Koeke rnoer et all 
I200I . We used a square interpolation kernel and set the 
output pixel scale at 0.05" pix^^. 

2.2.2. WFPC2 

Following the failure of ACS in the middle of Cycle 15 
observing, the Wide Field Camera CCD 3 of WFPC2 was 
used to image the remainder of the DOG population. For 
these observations, single-orbit data through the F606W 
filter were used to take advantage of WFPC2's superior 
sensitivity at this wavelength compared to other WFPC2 
filters. We used a four point dither pattern (WFPC2- 
BOX) with a point and line spacing of 0.559" and a 
pattern orientation of 26.6°. Total exposure time was 
Ril600 sec. The standard WFPC2 pipeline system was 
used to bias-subtract, dark-subtract, and flat-field the 
images (Mobasher et al., 2002). MultiDrizzle was then 
used to correct for geometric distortions, perform sky- 
subtraction, image registra tion, cosmic ray rejectio n and 
final drizzle combination ( Koekemoer et all l2002l ) . We 
used a square interpolation kernel and output pixel scale 
of 0.045" pix~^, leading to a per-pixel exposure time of 
Ri340 sec. Due to the irregular performance of WF4 and 
the PC CCDs, we have restricted our analysis to the 
WF2 and WF3 CCDs. 

2.2.3. NICMOS 

Single orbit data of the DOGs were acquired with 
NIC2 and the F160W filter. We used a two-point dither 
pattern (NIC-SPIRAL-DITH) with a point spacing of 
0.637". Total exposure time was w2600 s. To reduce the 
data, we followed the standard reduction process outlined 
in the NICMOS data handbook (McLaughlin & Wik- 
land 2007). We used the IRAF routine nicpipe to pre- 
process the data, followed by the biaseq task to correct 
for non-linear bias drifts and spatial bias jumps. We then 
used nicpipe a second time to do flat-fielding and initial 
cosmic-ray removal. The IRAF task pedsky was used to 
fit for the sky level and the quadrant-dependent residual 
bias. Significant residual background variation remained 
after this standard reduction process. To minimize these 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., u nder cooperative agree ment with 
the National Science Foundation. |http:/ /iraf.noao.edu/| 
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TABLE 1 
Observations 



Source Name 



RA (J2000) DEC (J2000) 



Optical Exposures 
Instrument/Filter UT Date 



Infrared Exposures 
Instrument/Filter UT Date IDS 



SST24 


J142538, 


.2+351855 


14: 


:25: 


:38, 


.155 


+35: 


:18: 


:56. 


.19 


2, 


.26=" 


WFPC2/F606W 


2007- 


-04- 


■29 


NIC2/F160W 


2007- 


■06- 


■16 


19 


SST24 


J142622, 


.0+345249 


14: 


:26: 


:22, 


.032 


+34: 


:52: 


:49. 


.69 


2, 


.00= 


ACS/F814W 


2006- 


■11- 


■25 


NIC2/F160W 


2007- 


■03- 


■13 


9 


SST24 


J142626, 


.4+344731 


14: 


:26: 


:26, 


.538 


+34: 


:47: 


:31. 


.53 


2, 


.13=" 


WFPC2/F606W 


2007- 


■12- 


■31 


NIC2/F160W 


2007- 


■01- 


■04 


16 


SST24 


J142644, 


.3+333051 


14: 


:26: 


:44, 


.321 


+33: 


:30: 


:52. 


.20 


3 


.312^ 


WFPC2/F606W 


2007- 


■04- 


■10 


NIC2/F160W 


2007- 


■02- 


■25 


31 


SST24 


J142645, 


.7+351901 


14: 


:26: 


:45, 


.701 


+35: 


:19: 


:01. 


.17 


1, 


.75=^ 


WFPC2/F606W 


2007- 


■04- 


■24 


NIC2/F160W 


2007- 


■05- 


■05 


2 


SST24 


J142648, 


.9+332927 


14: 


:26: 


:48, 


.970 


+33: 


:29: 


:27. 


.56 


2, 


.OO'^ 


ACS/F814W 


2007- 


■01- 


■17 


NIC2/F160W 


2006- 


■12- 


■19 


10 


SST24 


J142653, 


.2+330220 


14: 


:26: 


:53, 


.285 


+33: 


:02: 


:21. 


.37 


1, 


.86=^ 


ACS/F814W 


2006- 


■12- 


■29 


NIC2/F160W 


2007- 


■03- 


■01 


6 


SST24 


J142804, 


.1+332135 


14: 


:28: 


:04, 


.133 


+33: 


:21: 


:34. 


.97 


2, 


.34=^ 


WFPC2/F606W 


2007- 


■04- 


■17 


NIC2/F160W 


2007- 


■05- 


■01 


20 


SST24 


J142924, 


.8+353320 


14: 


:29: 


:24, 


.811 


+35: 


:33: 


:21. 


.30 


2, 


.73=^ 


WFPC2/F606W 


2007- 


■03- 


■18 


NIC2/F160W 


2007- 


■02- 


■13 


30 


SST24 


J142958, 


.3+322615 


14: 


:29: 


:58, 


.354 


+32: 


:26: 


:15. 


.17 


2, 


.64=^ 


WFPC2/F606W 


2007- 


■03- 


■14 


NIC2/F160W 


2007- 


■06- 


■22 


29 


SST24 


J143001, 


.9+334538 


14: 


:30: 


:01, 


.910 


+33: 


:45: 


:38. 


.54 


2, 


.46=^ 


WFPC2/F606W 


2007- 


■04- 


■12 


NIC2/F160W 


2007- 


■04- 


■28 


22 


SST24 


J143025, 


.7+342957 


14: 


:30: 


:25, 


.764 


+34: 


:29: 


:57. 


.29 


2, 


.545° 


WFPC2/F606W 


2007- 


■04- 


■26 


NIC2/F160W 


2007- 


■04- 


■13 


26 


SST24 


J143102, 


.2+325152 


14: 


31: 


:02, 


.220 


+32: 


:51: 


:52. 


.10 


2, 


.00'' 


WFPC2/F606W 


2007- 


■04- 


■26 


NIC2/F160W 


2008- 


■01- 


■07 


11 


SST24 


J143109, 


.7+342802 


14: 


31: 


:09, 


.823 


+34: 


:28: 


:02. 


.34 


2, 


.lO'^ 


WFPC2/F606W 


2007- 


■04- 


■12 


NIC2/F160W 


2007- 


■01- 


■04 


13 


SST24 


J143135, 


.2+325456 


14: 


:31: 


:35, 


.309 


+32: 


:54: 


:56. 


.84 


1, 


.48= 


WFPC2/F606W 


2007- 


■03- 


■21 


NIC2/F160W 


2007- 


■06- 


■17 


1 


SST24 


J143225, 


.3+334716 


14: 


:32: 


:25, 


.433 


+33: 


:47: 


:16. 


.67 


2, 


.00= 


ACS/F814W 


2006- 


■12- 


■07 


NIC2/F160W 


2007- 


■06- 


■17 


12 


SST24 


J143242, 


.5+342232 


14: 


32: 


:42, 


.569 


+34: 


:22: 


:32. 


.23 


2, 


.16=^ 


WFPC2/F606W 


2007- 


■04- 


■19 


NIC2/F160W 


2007- 


■12- 


■07 


18 


SST24 


J143251, 


.8+333536 


14: 


32: 


:51, 


.873 


+33: 


:35: 


:35. 


.89 


1, 


.78=^ 


WFPC2/F606W 


2007- 


■03- 


■20 


NIC2/F160W 


2008- 


■01- 


■14 


3 


SST24 


J143312, 


.7+342011 


14: 


33: 


:12, 


.734 


+34: 


:20: 


:11. 


.10 


2, 


.119^* 


WFPC2/F606W 


2007- 


■04- 


■20 


NIC2/F160W 


2007- 


■06- 


■19 


15 


SST24 


J143325, 


.8+333736 


14: 


:33: 


:25, 


.884 


+33: 


:37: 


:36. 


.90 


1, 


.90= 


WFPC2/F606W 


2007- 


■05- 


■01 


NIC2/F160W 


2007- 


■04- 


■20 


7 


SST24 


J143358, 


.0+332607 


14: 


:33: 


:58, 


.077 


+33: 


:26: 


:07. 


.46 


2, 


.414^ 


ACS/F814W 


2006- 


■12- 


■10 


NIC2/F160W 


2008- 


■01- 


■24 


21 


SST24 


J143447, 


.7+330230 


14: 


:34: 


:47, 


.762 


+33: 


:02: 


:30. 


.46 


1, 


.78=^ 


WFPC2/F606W 


2007- 


■03- 


■19 


NIC2/F160W 


2006- 


■12- 


■23 


4 


SST24 


J143504, 


.1+354743 


14: 


:35: 


:04, 


.166 


+35: 


:47: 


:43. 


.79 


2, 


.13=" 


WFPC2/F606W 


2007- 


■04- 


■26 


NIC2/F160W 


2007- 


■01- 


■03 


17 


SST24 


J143508, 


.4+334739 


14: 


:35: 


:08, 


.518 


+33: 


:47: 


:39. 


.44 


2, 


.10= 


WFPC2/F606W 


2007- 


■04- 


■12 


NIC2/F160W 


2007- 


■04- 


■22 


14 


SST24 


J143520, 


,7+340418 


14: 


:35: 


:20, 


.801 


+34: 


:04: 


:18. 


.30 


1, 


.79=^ 
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residuals, we constructed a normalized, object-masked 
median sky image based on all of our NIC2 science 
frames. This sky image was then scaled by a constant 
factor and subtracted from each science image. The scal- 
ing factor was computed by minimizing the residual of 
the difference between the masked science image and the 
scaled sky image. Mosaicing of the dithered exposures 
was performed using calnicb in IRAF, resulting in a 
pixel scale of 0.075" pix~^. 



Sky Survey frames, see lEisenhardt etldl l2004f ). Using 
the properly aligned, multi-wavelength dataset, identify- 
ing the proper counterpart to the MIPS source is rela- 
tively straightforward, since inspection of the four IRAC 
channels reveals a single source associated with the 24/im 
emission for all but one source (this source is undetected 
in all four IRAC channels). The absolute uncertainty in 
the centroid of the IRAC S.G/um emission ranges from 
0'.'3-0'.'5. 



2.3. Astrometry 

Each ACS/WFPC2 and NICMOS image is aligned to 
the reference frame of the NDWFS, which itself is tied 
to the USNO A- 2 catalog. We first run S ource Extrac- 
tor (SExtractor, 'Berti n fc ArnoutsI |1996[ ) on a cutout 
of the /-band NDWFS corresponding to the appropri- 
ate ACS/WFPC2 Field Of View (FOV) to generate a 
list of comparison objects. The IRAF task wcsctran is 
used to convert this list into pixel coordinates on the 
ACS/WFPC2 image. Another IRAF task, imcentroid, 
is used to improve the accuracy of the pixel coordinates. 
Finally, the IRAF task ccmap applies a first order fit 
to correct the zero point of the astrometry and update 
the appropriate WCS information in the header of the 
ACS/WFPC2 image. This ahgned ACS/WFPC2 image 
is then used as the reference frame for correcting the as- 
trometry of the NICMOS image and the IRAC images 
(since the IRAC images of the Bootes Field are not tied 
to the USNO A-2 catalog, but instead to the 2/xm AU- 



2.4. Photometry 

We perform 2" diameter aperture photometry on each 
DOG in both the rest-optical and rest-UV, choosing the 
center of the aperture to be located at the peak flux 
pixel in the NICMOS images. We remove foreground and 
background objects using SExtractor (see Section |4. 2. 2p 
and calculate the sky level using an annulus with an inner 
diameter of 2" and a width of 1". We found that in some 
cases (particularly those NICMOS images where signifi- 
cant residual non-linearities remained), the flux density 
radial profile did not flatten at large radii. When this 
occurred, we determined the appropriate sky value by 
trial-and-error. We computed the background level and 
photometric uncertainty by measuring the sigma-clipped 
mean and RMS of fluxes measured in N 2" diameter 
apertures, where « 10 and TV « 50 for the NICMOS 
and ACS/WFPC2 images, respectively. 

We compute 4" diameter aperture photometry in the 
NDWFS Bw^ R, and / images centered on the IRAC 
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3.6 /im centroid of emission. Sky background levels were 
computed in a 3" wide annulus with an inner diameter of 
4". Limiting magnitudes were determined by measuring 
the flux within a 4" aperture at several sourceless loca- 
tions near the DOG and computing the rms variation of 
the flux values. 

We verified the accuracy of our ACS and WFPC2 pho- 
tometric zeropoints by comparing well-detected sources 
common to both our HST and NDWFS imaging. For our 
ACS/F814W observations, we compared to the NDWFS 
/-band imaging and found negligible offsets (-0.03 ± O.fO 
magnitudes). For our WFPC2/F606W observations, we 
compared to the NDWFS i?-band (after correcting for 
color terms due to the dissimilarity of the R and F606W 
filter bandpasses) and again found negligible offsets (0.05 
± 0.15 magnitudes). 

2.5. Images of DOGs 

Figure [3] shows 2" x 2" cutout images of the DOGs 
in order of increasing rcdshift. Each cutout is centered 
roughly on the centroid of emission as seen in the NIC- 
MOS image. A red plus sign shows the centroid of IRAC 
3.6/im emission and is sized to represent the l-cr uncer- 
tainty in the position, which includes independent contri- 
butions from the centroiding error on the 3.6/im emission 
(«0'.'2-0'.'4, depending on S/N), the relative astrometric 
calibration uncertainty within the 3.6/im map («0'.'2), 
and the uncertainty in tying the 3.6/im map to the HST 
images («0'.'l). The la rms offset between IRAC and 
NICMOS centroids of the sample is 0'.'2. In most cases, 
the offset in centroids is negligible, but those cases where 
it is not are associated with faint 3.6/im emission (when 
the absolute astrometric uncertainty may be as large as 
0'.'5). This suggests there is no significant offset between 
the near-IR and mid-IR centroids, although we note that 
we cannot rule out offsets at the < 1 kpc scale. 

The DOGs exhibit a wide range of morphologies, 
with most being well-resolved. Only one object (SST24 
J142644. 3-1-333051) shows strong Airy rings and is clearly 
an unresolved point source. Here we give a brief qualita- 
tive desription of the morphology of each object. 

(1) SST24 J143135. 2+325456: F606W: Weak de- 
tection. F160W: Large-scale emission with a faint tail 
extending northeast. 

(2) SST24 J142645. 7+351901: F606W: No signif- 
icant detection. F160W: Two compact, resolved compo- 
nents separated by w0'.'5. 

(3) SST24 J143251. 8+333536: F606W: Faint 
emission slightly NE of NICMOS centroid; possible sec- 
ond source --0'.'5 northwest (NW) of NICMOS centroid. 
F160W: Extended object; possible point source contam- 
ination. 

(4) SST24 J143447.7+330230: F606W: No signif- 
icant detection. F160W: Irregular, diffuse object. 

(5) SST24 J143520.7+340418: F606W: Compact, 
resolved object. F160W: Very compact object, but no 
evidence for Airy rings. 

(6) SST24 J142653. 2+330220: F814W: Large 
scale, irregular, and diffuse emission. F160W: Large- 
scale, irregular, and diffuse, but with bright compact nu- 
clear component. 

(7) SST24 J143325.8+333736: F606W: Compact, 
resolved object. F160W: Bright compact and extended 
components. 



(8) SST24 J143644.2+350627: F606W: Compact, 
resolved object. F160W: Extended object; possible point 
source contamination . This object has a counterpart in 
the XBootes catalog (jBrand et al.ll2006D . 

(9) SST24 J142622. 0+345249: F814W: Four com- 
pact, resolved clumps spread in a 'T' shape with no vis- 
ible central component. F160W: Similar irregular 'T' 
shape, but components are not as distinct. NE compo- 
nent is bluer than other components. 

(10) SST24 J142648. 9+332927: F814W: Com- 
pact, resolved object. Chain of sources extends towards 
southwest. One of these sources is within 1" of the DOG 
and is included in the photometric and morphological 
measurements, since there is no clear evidence to sug- 
gest it is not associated with the system. F160W: Similar 
compact, resolved object; possible point source contam- 
ination. 

(11) SST24 J143102. 2+325152: F606W: No signif- 
icant detection. F160W: No significant detection. This 
object is also not detected in any of the IRAC images, 
but does have a counterpart in the XBootes catalog 
(jBrand et al.ll2006D . 

(12) SST24 J143225. 3+334716: F814W: Com- 
pact, irregular object. F160W: Compact, resolved ob- 
ject; possible point source contamination. 

(13) SST24 J143109.7+342802: F606W: No sig- 
nificant detection. F160W: Irregular, diffuse object. 

(14) SST24 J143508.4+334739: F606W: Com- 
pact, resolved, and irregular central component with tail 
of emission to southeast (SE). F160W: Very similar, but 
central component is stronger relative to tail. 

(15) SST24 J143312. 7+342011: F606W: Four 
compact components in a semi-circle offset from the cen- 
troid of NICMOS emission. F160W: Extended object; 
possible point source contamination. 

(16) SST24 J142626.4+344731: F606W: No sig- 
nificant detection. F160W: Irregular, diffuse object. 

(17) SST24 J143504.1+354743: F606W: Faint 
source barely detected. F160W: Irregular, extended ob- 
ject; possible point source contamination. 

(18) SST24 J143242. 5+342232: F606W: No sig- 
nificant detection. F160W: Faint, compact compo- 
nent near 3.6/im centroid with emission leading to sec- 
ond, brighter peak '-~^0'.'5 to southwest (SW). Possible 
multiple-component system. 

(19) SST24 J142538.0+332607: F606W: No sig- 
nificant detection. F160W: Irregular objected elongated 
in NW-SE direction. 

(20) SST24 J142804.1+332135: F606W: No sig- 
nificant detection. F160W: Faint, irregular object. 

(21) SST24 J143358.0+332607: F606W: Faint, ir- 
regular object. F160W: Extended object with possible 
point source contamination. 

(22) SST24 J143001.9+334538: F606W: No sig- 
nificant detection. F160W: Faint, irregular object. 

(23) SST24 J143545. 1+342831: F814W: Faint, 
compact irregular objects; possibly two component sys- 
tem. F160W: Extended emission; possible point source 
contamination. 

(24) SST24 J143725. 1+341502: F814W: Very 
faint, low surface brightness feature extending to east. 
F160W: Diffuse emission with compact object; faint Airy 
ring present. 

(25) SST24 J143808. 3+341016: F814W: Faint, 
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Fig. 3. — Cutouts of the 31 DOGs observed by HST, shown with a linear stretch. Columns 1 and 3 are the rest-UV images from either 
ACS F814W or WFPC2 F606W. Columns 2 and 4 are the rest-optical images from NIC2 F160W. Each cutout is 2" on a side and is 
oriented north up and east left. The objects are arranged in order of increasing redshift, and the redshift is printed in the lower right corner 
of each NICMOS image. A red cross denotes the position and l-cr uncertainty in the centroid of the IRAC 3.6/im emission. In NICMOS 
images where the S/N per pixel is greater than 2, white contours outline the brightest 20% pixels, and black contours show the outline of 
the segmentation map used in measuring the non-parametric morphologies. 



compact, resolved components offset from centroid of 
N1C2 emission. F160W: Compact, central component 
with extension to SE overlapping ACS emission centroid. 

(26) SST24 J143025. 7+342957: F606W: Com- 
pact, resolved object. F160W: Compact, resolved object; 
possible point source contamination. 

(27) SST24 J143523. 9+330706: F814W: Com- 
pact, resolved object with tail of emission to SE; possible 
point source contamination. F160W: Compact, resolved 
object; possible point source contamination. 

(28) SST24 J143539. 3+334159: F606W: Possible 
faint diffuse emission N of NIC2 centroid. F160W: Com- 
pact, resolved object; possible point source contamina- 



tion; possible tail of emission towards N. 

(29) SST24 J142958.3+322615: F606W: Compact 
resolved object; possible point source contamination. 
F160W: Extended object with bright nuclear source; pos- 
sible point source contamination. 

(30) SST24 J142924.8+353320: F606W: No signif- 
icant detection. F160W: Very compact, irregular, faint 
object near IRAC 3.6/im centroid. Larger, brighter ob- 
ject --O'.'S to north. 

(31) SST24 J142644.3+333051: F606W: Weak de- 
tection. F160W: Dominated by point source emission; 
clear Airy ring. This so urce has two X-ray counterparts 
in the XIBootes catalog ([Brand et al.li200a ). 



HST PL DOG Morphologies 



(9) SST24 J142622.0+345249 



(10) SST24 J142648.9+332927 




NIC2/F160W 




(11) SST24 J143 102 2+325152 
"wft^CZ/i^ieW ? n I^C2/F160W 



(12) SST24 J143225.3+334716 




^STfsMwvF ^y 
■. ■. - iJ- ■ % 



NIC2/F160W 




2.!)0 



(13) SST24 J143 10 9.7+342802 



(14) SST24 J143508.4+334739 



(15) SST24 J1 433 12.7+342011 
WpPeE/F606W 

S-1. 




NIC2/F160W 




2X19 





NIC2/F1;60W 














L 




^^^.2.10 





















(16) SST24 J142626.4+344731 




Fig. 3. — Continued. 

3. METHODOLOGY: MORPHOLOGICAL 
ANALYSIS 

We undertook three different, complementary ap- 
proaches to analyzing the morphology of the DOGs in 
our sample: a visual classification experiment, multi- 
component GALFIT modeling, and non-parametric 
quantification. In this section, we describe the details of 
our methodology. The results are described in section [H 

3.1. Visual Classification 

We first undertook a visual classification of the 
DOGs by conducting the following experiment: for each 
ACS/WFPC2 image, we generate a 5"x5" cutout image 
of both the DOG and 14 other randomly selected galax- 
ies in the same FOV with the same magnitude range as 
our DOG sample. Eight of the coauthors then classified 



all 15 galaxies in each FOV (the DOG was not identi- 
fied), placing them into one of the following bins: Ellip- 
tical/Compact (E/C), Disk, Irregular/Multi-component, 
Irregular/Diffuse, or Too Faint to Tell. In practice, since 
these galaxies are selected to be faint in the optical and 
generally have low S/N, we group together the E/C and 
Disk categories into a "Regular" bin and the two irregular 
categories into an "Irregular" bin. This results in a total 
of 3600 independent classifications, of which 240 pertain 
to the DOGs. In an effort to explore the robustness of 
our results, we fiag and remove from the sample all ob- 
jects where fewer than six classifers were in agreement. 
This has the effect of reducing the fraction of "Too Faint 
To Tell" responses, but the ratio of the Regular to Irreg- 
ular classifications changes by less than 15%. A similar 
experiment was done on the NICMOS images, without 
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Fig. 3. — Continued. 



the control sample, since the NIC2 FOV is so small. 

Interpretation of the results of our visual classifica- 
tion analysis is hampered by low S/N (in the case of 
the ACS/WFPC2 images) or the lack of a control sam- 
ple (in the case of the NICMOS images), so we forego 
a detailed analysis and instead present the mode of the 
classifications for each DOG along with an indication of 
whether the eight coauthors were in general agreement 
in our morpholical results table in section 14.2 21 This 
is useful as a qualititative assessment of the morphol- 
ogy for comparison with the more quantitative methods 
discussed below. 

3.2. GALFIT Modeling 

In many of the NICMOS images, there is a com- 
pact component that is not seen in the corresponding 



ACS/WFPC2 image, implying there is significant obscu- 
ration in the central region of many of the DOGs. In 
this section, we describe the method we use to explore 
the degree to which each DOG is dominated by a cen- 
tral, unre solved componen t. Our tool in this effort is 
GALFIT fPcn g et al.[|2002D . which uses a 2-dimensional 
minimization to search the parameter space of a set 
of predefined functions and identify the parameters that 
best describe the observed 2-D profile. 

Because the DOGs are small and have low S/N com- 
pared to more typical applications of GALFIT, we re- 
strict the size of the fitting region to be 41x41 pixels 
(corresponding to an angular and physical size of 3" and 
~24 kpc, respectively) and include the minimum neces- 
sary components in our model. For a variety of reasons, 
we expect AGN to be important contributors to the emit- 
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ted radiation from these sources. Therefore, we model 
the observed emission with three components which are 
described by a total of 10 free parameters. The number 
of degrees of freedom, iVnoF , is calculated as the number 
of pixels in the image being modeled minus the number of 
free parameters. This implies that the maximum A^dof 
is 1671. Those cases where A'dof < 1671 are associated 
with images where pixels were masked out because they 
were associated with residual instrumental noise and pre- 
vented convergence with GALFIT. We note that because 
NIC2 is a Nyquist-sampled array (0.075" pix~^ compared 
to 0.16" FWHM beam), the pixels in our image may not 
be completely independent. As a result, the values 
should be interpreted in a relative sense rather than an 
absolute one. 

The first element in our GALFIT model is a sky com- 



ponent whose amplitude is chosen to obtain flat radial 
profiles at large radii and is not allowed to vary. The sec- 
ond is an instrumental PSF generated from the TinyTim 
software (Krist and Hook 2004), which can simulate a 
PSF for NICMOS, WFPC2, and ACS. For the NICMOS 
and WFPC2 images, the DOG is positioned in nearly the 
same spot on the camera. In the case of WFPC2 this is 
pixel (400,400) of chip 3 and pixel (155, 164) for NIC- 
MOS. Meanwhile, a different region of the ACS camera 
is used for each DOG. Therefore, we generate a unique 
PSF at each position on the appropriate chip in which 
a DOG is observed. We use a red power-law spectrum 
{Fy oc v^^) as the object spectrum. The PSF is com- 
puted out to a size of 3.0", and for the WFPC2 PSF we 
oversample by a factor of 2 to match the pixel scale of 
the drizzled WFPC2 images. 
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The final component is a Sersic profile (jSersid [l968f) 
where the surface brightness scales with radius as 
exp[-K((r/i?cff)^/"-l)], where K is chosen such that half 
of the flux falls within i?cff- We attempted to place as 
few constraints as possible so as to optimize the measure- 
ment of the extended flux (i.e., non-point source compo- 
nent). However, in certain cases, the Sersic index had 
to be constrained to be positive to ensure convergence 
on a realistic solution. For the NICMOS images, we 
used the uncertainty image output by calnicb as the 
error image required by GALFIT to perform a proper 
minimization. The Tiny Tim NIC2 PSF is convolved 
with the Sersic profile prior to performing the mini- 
mization. The initial guesses of the magnitude, half-light 
radius, position angle, and ellipticity were determined 
from the output values from SExtractor. Varying the ini- 
tial guesses within reasonable values (e.g., plus or minus 
two pixels for the half-light radius) yielded no significant 
change in the best-fit model parameters. We used the 
NICMOS centroid as the initial guess for the (x,y) posi- 
tion of both the PSF and extended components, but in a 
few cases these guesses had to be modified by 1-2 pixels 
in order to result in convergence. 

We note that we tested two-component models as well 
(single component Sersic profile plus sky background) 
and found larger reduced values, especially when 
the point source fraction in our three-component model 
was large (see further discussion in section 14.2. ip . In 
cases where the point source fraction was small, the two- 
component model had similar parameter values as the 
three-component model, as we would expect. 

It is important to note here that NIC2 cannot spatially 
resolve objects smaller than 1.3 kpc at z ~ 2. This limit 
is large enough to encompass a compact stellar bulge as 
well as an active galactic nucleus, implying that we can- 
not, from these data alone, distinguish between these two 
possibilities as to the nature of the central, unresolved 
component. 

After the best-fit parameters are found in the NICMOS 
image, we run GALFIT with a simplified model on the 
DOGs in the ACS/WFPC2 images. The primary simpli- 
fication is to fix the position of the PSF based on the best- 
fit NICMOS PSF location (allowing up to 2 pixel wiggle 
room to account for astrometric uncertainties, which can 
be as large as O'.'l). In many cases, GALFIT required an 
upper limit to be placed on the magnitude of the PSF 
component in order to reach convergence. We choose to 
use the magnitude of a point source detected at the 2-a 
level for this upper limit. We note that our Sersic profile 
model for the extended DOG flux is not representative 
of the rest-UV morphology of many of the DOGs (i.e., 
the reduced x^ values are large), but it does adequately 
recover their total flux. 

3.3. Non-parametric Classification 

The Gini coefficient (G) and A/20 parameter are known 
to be reliable tools for the chara cterization of faint-object 
morphologies (|Lotz et al.l l2004f l . G was originally cre- 
ated to mea sure how even ly the wealth in a society is 
distributed (iGlasser 1 119621). Recently, [Abraham et al.l 
dJOOS) and lLotz et all (|2004[ ) applied this method to aid 
in the classification of galaxies, with G defined such that 
low (high) values imply an equal (un-equal) distribution 
of flux. A/20 is the logarithm of the second-order mo- 



ment of the brightest 20% of the gala xy's flux, norma l- 
ized by the total second-order moment (jLotz et al.ll2004l ). 
This means that higher values of M20 imply multiple 
bright clumps offset from the second-order moment cen- 
ter. Lower values, on the other hand, suggest a system 
dominated by a central component. 

Prior to computing G or M20, we first generate a cat- 
alog of objects using SExtractor. We use a detection 
threshold of 1.5a (corresponding to 24.5 mag arcsec"^) 
and a minimum detection area of 15 pixels. The center 
of the image as well as the ellipticity and position angle 
computed by SExtractor are used as inputs for comput- 
ing morphological measures. In addition, we use catalog 
sources selected to have magnitudes within the range of 
all 24 DOGs analyzed in this paper as a "field" galaxy 
sample for comparison to DOGs. 

Much of the methodology in this section relies on mor- 
ph ology code writte n by J. Lotz and described in detail 
in iLotz et al.l ()2004[ ). Here, we summarize the relevant 
information. Postage stamps of each object in the SEx- 
tractor catalog (and the associated segmentation map) 
are created with foreground/background objects masked 
out. Using a small region of the cutout devoid of sources, 
a sky value is computed and subtracted from the postage 
stamp. Next, we determine which pixels in each postage 
stamp belong to the galaxy and which do not. Since the 
isophotal-based segmentation map produced by SExtrac- 
tor is subject to the effects of surface brightness dimming 
at high redshift, we use a segmentation map based on the 
mean surface brightness at the Petrosian radius /i(/?p). 
Pixels with surface brightness above /i(i?p) are assigned 
to the galaxy while those below it are not. We define Rp 
as the radius at which the ratio of the surface brightness 
at Rp to the mean surface brightness within Rp is equal 
to 0.2. 

Using the new segmentation map, we recompute the 
galaxy's center by minimizing the total second-order mo- 
ment of the flux. A new value of Rp is then computed 
and a revised segmentation map is used to calculate G 
and M20. Finally the morphology code produces an av- 
erage S/N per pixel value using the pi xels in the revise d 
segmentation map (Eqs. 1 through 5 in lLotz et al.ll200^ . 

One of the most common methods of characterizing 
galaxy morphologies in the literature is to rn easure the 
concentration index C ('Abraham et al."1 994[ ). the rota- 
tional asymmetry A (Schadc ct al. 1995), and the resid- 
ual dumpiness, S (jConselicel 120031 ) . Given sufficiently 
high S/N and spatial resolution, the GAS system has 
had demonstrated success in measuring mor phological 
param eters and identifyi ng mergers at low (IConselicd 
I2003D and high redshift (|Conselice et"alll2008f l. Unfor- 
tunately, the objects in our sample do not meet simulta- 
neously the S/N and spatial resolution requirements to 
be reliably placed in CAS space. Because computation 
of A and S involves differencing two images, the neces- 
sary per-pixel S/N to measure these parameters reliably 
is twice as high as those that do not involve subtracting 
images. We find per-pixel S/N values ranging from ^2-5 
for the DOGs, whereas reliable m easurements of A and S 
require S/N>5 (|Lotz et al.ll2004D . In principle, the data 
are of sufficient quality to measure G (see Tab.E]), but in 
practice we find that the inherent assumption of circular 
symmetry does not apply well to the DOGs, making the 
interpretation of C values difficult. 
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4. RESULTS 

4.1. Photometry 

In Figure [H we show the color-magnitude diagram for 
DOGs and a sample of galaxies in the HDF whose photo- 
metric redshifts are comparable to DOGs (1.5 < Zphot < 
2.5). For DOGs where the measured flux is below the 2(t 
detection limit, we use an open plotting symbol and an 
upward pointing arrow. DOGs range in H-hand magni- 
tude from 21.93 to 25.1 AB mags. In both V - H and 
I — H, DOGs are redder than a typical high-z galaxy 
by 0.2-3 AB mags. In partic ular, the LBGs from the 
HDF-N (|Papovich et al.ir200l[ ) are comparably bright in 
H, but fainter in V by ~2 AB mags than DOGs. There 
is substantial overlap between the colors of DOGs and 
DRGs, suggesting that we might expect to see similari- 
ties in the morphologies between these two populations. 
Table [2] summarizes the photometric information derived 
from the NDWFS and the HST imaging. 

4.2. Morphologies 
4.2.1. GALFIT Results 

The results of our GALFIT analysis for the extended 
component Sersic profile fit to the NICMOS images are 
shown in Table [31 along with l-cr uncertainties in the 
best-fit parameters. The Sersic indices (n) range from 
0.1 to 2.2 (median n — 0.9). For those objects where 
n < 1, we note that constraining n to be equal to 1 does 
not significantly alter the remaining fit parameters. As 
the Sersic index decreases, the radial profile flattens more 
rapidly within r < R^g and the intensity drops more 
steeply beyond r > i?eff- For reference, n = 0.5 corre- 
sponds to a Gaussian profile, n = 1 corresponds to an 
exponential prof ile, and n = 4 cor responds to a de Vau- 
couleurs profile (Peng et al."2002^. The ratio of the mi- 
nor to major axis ranges from 0.20 to 0.88 with a median 
value of 0.53. In comparison, simulated merger remnants 
tend to have a luminous component in the shape of an 
oblat e spheroid, with axis ratios of 1:1:0.5 (jNovak et al.l 
|2006[) . The projected axial ratio should thus vary be- 
tween 0.5 and 1.0. Our observed median value of wO.5 
suggests that DOGs have more disk-like profiles than the 
simulated merger remnants. This may be due to a non- 
merger origin for DOGs, or it may be an indication that 
DOGs have not progressed to the merger remnant stage. 

It is possible for a degeneracy to arise in the fitting 
parameters, in the sense that a high-n, low point source 
fraction model may be comparable to a low-n, high point 
source fraction model. We have tested this by running 
GALFIT with the point source fraction set to zero (i.e., 
removing the PSF component). The resulting n values 
range from 0.1 to 5.2 with a median of 2.0, which is 
still below the value of 4 that is typical of early-type 
galaxies. A total of 7 DOGs have n > 3 using this zero- 
point-source model. The best-fit i?off values change by 
less than one pixel, with an offset of -0.2±0.8 pixels. In 
general, as mentioned in section [3?2l the removal of the 
PSF component leads to larger reduced-x^ values (see 
columns 5 and 6 in Tab. [3]). However, we note that only 
one case (SST24 J142644.3-K333051) is associated with a 
>0.08 decrease in xt after adding a non-zero PSF com- 
ponent. This suggests that the PSF component in most 
of the DOGs in this sample is not dominant at rest-frame 
optical wavelengths. 



We find the effective radius, Rcs, ranges from 1.1 to 
5.9 kpc with a median value of 2.5 kpc. In the left panel 
of Figure O we show the distribution of i?eff values for 
DOGs (dark shaded histogram). Distant Blue Galax- 
ies (diagonal blue hatched; DBGs, i .e., galaxies wit h 
Zphot > 2 that satisfy J - Ks < 2.3; iToft et al.ll2007fl 
and DRGs (opposite diagonal red hatched: IZirm et all 
120071: iToft et al.l 120071 ). Two-sided K-S tests show that 
DOGs are dissimilar from both populations, with a <7% 
and <4% chance of being drawn from the same par- 
ent distribution, respectively. Based on the nature of 
their UV-NIR SED, DRGs may be separated into those 
that are actively forming stars (active DRGs or sDRGs) 
and those that a r e not (quiescent DRGs or qDRGs, see 
IZirm et al.l l2007t iToft et al.|[2007f ). The right panel of 
Figure [5] shows the DOG i?eff distribution in comparison 
to active DRGs (diagonal blue hatched) and quiescent 
DRGs (opposite diagonal red hatched) . Quiescent DRGs 
have much smaller effective radii, while active DRGs are 
much closer to DOGs. Here the two-sided K-S test gives 
a 99.9994% and 34% chance of being drawn from differ- 
ent parent distributions, respectively, suggesting that the 
DOG and active DRG populations may overlap. 

In Table [H we give the V, I, and H magnitudes of 
the nuclear (PSF) component and the extended (galaxy) 
component, as well as the fraction of light contributed 
by a point source (including 1-a uncertainties). When 
the nuclear component is not detected, we quote the 3- 
(7 limit on the point source fraction. The magnitude of 
the PSF component is measured using the same aper- 
ture photometry method described in Section 12.41 with 
the exception that the sky background is assumed to be 
zero. For the extended component, we subtract the PSF 
component from the science image and compute the pho- 
tometry in the usual way on the residual image. The 
fraction of light due to an unresolved component in the 
rest-optical ranges from 0.04 to 0.78, and the median 
is 0.12. In the rest-UV, however, this fraction is signi- 
cantly smaller, with only one object having a detected 
fraction. This object, SST24 J143523.9-H330706, stands 
out as unique by virtue of having a greater point source 
fraction in the rest-UV than in the rest-optical. This be- 
havior is unique within our sample (but is expected when 
the AGN is viewed without obscuration) and is also re- 
flected in the non-parametric measures of its morphology 
(see Sect. 14 .272] for more detail). 

Figure [6] shows the V — H and I — H colors of the 
nuclear, extended, and full galaxy components as a func- 
tion of H, in AB magnitudes. High-z galaxies and DRGs 
in the HDF-N and HDF-S are also shown. The full 
galaxy and extended components have similar colors and 
H magnitudes, consistent with the nuclear component 
not dominating the flux. This is why even the extended 
components of DOGs are redder in both V — H and I — H 
compared to Lyman break galaxies (LBGs) in the HDF, 
and suggests that one cannot create a DOG simply by 
adding an obscured AGN to a star-forming galaxy like 
an LBG. DRGs show greater overlap with the colors of 
DOGs but few are as bright in H as the DOGs in our 
sample. 

4.2.2. Non-parametric Classification Results 

Non-parametric methods of characterizing galaxy mor- 
phology are known to require high S/N imaging to yield 
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Fig. 4. — Color-magnitude diagram for DOGs. Left: V — H vs. H for DOGs observed by WFPC2 (filled black squares show 
detections, open black squares show lower limits). Galaxies spanning the redshift range 1.5 < z < 2.5 in the HDF-N (Papovich, personal 
communication) and HDF-S fjLabbe et al.ii200S^ are shown with black dots. Bright LBGs from the HDF-N are shown with green diamonds 
(jPapovich ct al. 2001). DRGs in the HDF-S are represented by filled orange triangles. Right: I — H vs. H for DOGs observed by ACS. 
Symbols are the same as in left panel. 



TABLE 2 
Photometric Properties'^'' 



Source Name 



Bw 



R 



y {F606W) / (F814W) H (F160W) F24 (mjy) R ~ [24]^ 



SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 
SST24 



J142538. 
J142622. 
J142626. 
J142644. 
J142645. 
J142648. 
J142653. 
J142804. 
J142924. 
J142958. 
J143001. 
J143025. 
J143102. 
J143109. 
J143135. 
J143225. 
J143242. 
J143251. 
J143312. 
J143325. 
J143358. 
J143447. 
J143504. 
J143508. 
J143520. 
J143523. 
J143539. 
J143545. 
J143644. 
J143725. 
J143808. 



2+351855 
0+345249 
4+344731 
3+333051 
7+351901 
9+332927 
2+330220 
1+332135 
8+353320 
3+322615 
9+334538 
7+342957 
2+325152 
7+342802 
2+325456 
3+334716 
5+342232 
8+333536 
7+342011 
8+333736 
0+332607 
7+330230 
1+354743 
4+334739 
7+340418 
9+330706 
3+334159 
1+342831 
2+350627 
1+341502 
3+341016 



> 26.6 
24.5+0.1 

> 26.6 

> 26.5 

> 26.6 
25.1+0.2 

> 26.6 

> 26.4 

> 26.6 
25.6+0.1 

> 26.4 
24.6+0.1 

> 25.7 

> 26.4 
24.7+0.1 

> 26.9 

> 26.5 

> 26.5 
24.5+0.1 
25.6+0.3 

> 26.7 

> 26.6 

> 26.6 
24.6+0.1 
24.8+0.1 

> 26.8 

> 26.4 

> 26.9 
24.7+0.1 
25.4+0.2 
25.1+0.1 



> 25.9 
24.5+0.3 

> 25.4 
24.3+0.2 

> 25.8 

> 25.0 

> 26.1 

> 25.7 

> 25.4 

> 25.7 

> 25.8 
24.0+0.1 

> 25.2 

> 25.5 
23.9+0.1 

> 25.3 

> 25.3 

> 25.5 
24.2+0.2 
24.6+0.3 

> 25.9 

> 26.1 

> 25.8 
24.1+0.1 

> 25.1 
25.0+0.2 

> 25.5 

> 25.2 
24.3+0.1 

> 25.4 
24.4+0.3 



> 25.5 
24.0+0.3 

> 25.2 
24.3+0.2 
24.5+0.3 
24.1+0.1 
24.7+0.3 

> 25.3 

> 24.9 

> 25.4 

> 25.1 
23.9+0.1 

> 25.2 

> 25.2 
23.5+0.1 

> 25.2 

> 25.1 

> 25.1 
23.9+0.1 
23.9+0.2 

> 25.3 

> 25.2 

> 25.5 
23.8+0.1 
24.1+0.2 
24.9+0.3 
24.7+0.3 

> 25.0 
24.0+0.2 

> 25.2 
23.9+0.1 



> 26.0 

> 26.4 
25.9+0.3 

> 26.2 



> 26.6 

> 26.1 
25.5+0.3 

> 26.5 
24.21+0.07 

> 26.0 

> 26.3 
25.1+0.2 

> 26.7 

> 26.5 
24.7+0.1 
25.6+0.3 

> 26.0 

> 26.4 
23.87+0.05 

> 26.2 

25.1+0.2 
25.6+0.2 



24.18+0.06 



24.9+0.2 
25.0+0.2 



26.0+0.4 

> 25.9 

24.7+0.2 

> 26.7 

26.1+0.6 
24.7+0.2 



24.0+0.1 
23.6+0.1 
23.7+0.1 
21.93+0.02 
23.31+0.09 
23.3+0.1 
22.7+0.1 
25.1+0.5 
24.7+0.3 
23.26+0.09 

24.9+0.3 
22.29+0.03 
> 25.1 
23.6+0.1 
22.04+0.03 

22.8+0.1 
22.68+0.04 
22.20+0.03 
22.23+0.04 
21.52+0.03 
23.09+0.06 

23.1+0.1 
23.08+0.09 
22.69+0.07 

24.0+0.2 
22.93+0.06 

23.1+0.1 
22.59+0.05 
22.70+0.07 
22.70+0.08 
22.34+0.08 



0.85+0.05 
1.29+0.05 
1.17+0.04 
1.14+0.04 
1.14+0.05 
2.33+0.07 
0.88+0.05 
0.84+0.03 
1.04+0.05 
1.18+0.05 
3.84+0.06 
2.47+0.05 
1.19+0.05 
1.11+0.04 
1.51+0.05 
1.28+0.05 
0.91+0.04 
0.82+0.04 
1.76+0.04 
1.87+0.06 
1.07+0.04 
1.71+0.04 
1.26+0.05 
2.65+0.08 
1.53+0.06 
1.09+0.05 
2.67+0.06 
1.95+0.05 
2.34+0.05 
1.41+0.05 
1.71+0.05 



> 16.1 
15.2 

> 16.0 
14.9 

> 16.3 

> 16.3 

> 16.3 

> 15.9 

> 15.9 

> 16.2 

> 17.7 
15.4 

> 15.8 

> 16.0 
14.8 

> 16.0 

> 15.6 

> 15.7 
15.2 
15.6 

> 16.3 

> 17.1 

> 16.5 
15.6 

> 15.9 
15.5 

> 16.9 

> 16.3 
15.6 

> 16.2 
15.4 



b c 
magnitude lower limits represent 2(7 values. magnitudes given in AB system. R — [24] color in Vega system. 
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TABLE 3 
GALFIT Results 



RcS 

n Axial Ratio (kpc) xl xl'^ ^^dof 



SST24 J142538.2+351855 0.7±0.2 

SST24 J142622.0+345249 O.liO.l 

SST24 J142626.4+344731 0.8±0.2 

SST24 J142644.3+333051 5.6±3.4 

SST24 J142645.7+351901 O.liO.l 

SST24 J142648.9+332927 0.7±0.3 

SST24 J142653.2+330220 0.4±0.1 

SST24 J142804.1+332135 1.0±0.6 

SST24 J142924.8+353320 0.8±0.3 

SST24 J142958.3+322615 0.4±0.1 

SST24 J143001. 9+334538 1.0±1.3 

SST24 J143025.7+342957 l.OiO.l 

SST24 J143102.2+325152 — 

SST24 J143109.7+342802 1.2±0.4 

SST24 J143135.2+325456 0.5±0.1 

SST24 J143225.3+334716 0.9±0.2 

SST24 J143242.5+342232 2.0±0.4 

SST24 J143251. 8+333536 0.8±0.1 

SST24 J143312.7+342011 0.7±0.1 

SST24 J143325.8+333736 1.3±0.1 

SST24 J143358.0+332607 2.1+0.5 

SST24 J143447.7+330230 0.4+0.1 

SST24 J143504.1+354743 0.5+0.1 

SST24 J143508.4+334739 2.6+0.4 

SST24 J143520.7+340418 0.6+0.4 

SST24 J143523.9+330706 0.5+0.1 

SST24 J143539.3+334159 1.8+0.5 

SST24 J143545.1+342831 1.1+0.2 

SST24 J143644.2+350627 0.9+0.2 

SST24 J143725.1+341502 0.3+0.1 

SST24 J143808.3+341016 0.9+0.1 



0.64+0.09 


2.5+0.4 


0.35 


0.36 


1653 


0.55+0.04 


2.5+1.0 


0.92 


0.92 


1671 


0.59+0.07 


3.0+0.4 


1.35 


1.35 


1671 


0.71+0.07 


1.1+0.3 


1.16 


2.58 


1671 


0.41+0.03 


4.6+0.2 


1.08 


1.08 


1671 


0.70+0.05 


1.8+0.1 


1.08 


1.08 


1561 


0.85+0.03 


2.9+0.1 


0.95 


0.96 


1671 


0.19+0.06 


6.4+2.6 


0.95 


0.96 


1671 


0.35+0.06 


2.8+0.4 


1.04 


1.07 


1671 


0.78+0.03 


2.1+0.1 


0.36 


0.39 


1671 


0.6+0.1 


1.2+0.3 


0.49 


0.49 


1671 


0.52+0.02 


1.8+0.1 


0.90 


0.96 


1671 


0.36+0.05 


7.0+1.7 


1.15 


1.15 


1671 


0.53+0.01 


3.8+0.1 


1.54 


1.56 


1671 


0.55+0.03 


1.9+0.1 


0.93 


0.94 


1671 


0.50+0.06 


4.0+0.8 


0.57 


0.58 


1601 


0.50+0.02 


3.3+0.1 


0.50 


0.52 


1671 


0.51+0.01 


4.0+0.1 


0.63 


0.66 


1671 


0.70+0.01 


3.1+0.1 


0.52 


0.60 


1671 


0.84+0.07 


1.4+0.1 


0.49 


0.48 


1671 


0.88+0.04 


1.9+0.1 


0.88 


0.89 


1671 


0.43+0.03 


4.6+0.3 


0.91 


0.93 


1671 


0.40+0.03 


2.2+0.2 


1.06 


1.10 


1671 


0.5+0.1 


3.4+0.7 


1.08 


1.10 


1671 


0.46+0.02 


1.5+0.1 


1.00 


1.05 


1671 


0.82+0.09 


3.7+0.7 


0.91 


0.92 


1671 


0.49+0.03 


2.2+0.1 


1.04 


1.06 


1671 


0.70+0.05 


1.9+0.1 


1.22 


1.23 


1670 


0.50+0.03 


3.4+0.2 


1.09 


1.14 


1424 


0.81+0.04 


3.0+0.2 


1.37 


1.38 


1671 



^ Asssuming finite PSF contribution Asssuming zero PSF contribution 
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Fig. 5. — Left: Distribution of effective radius, Rcfti for an unconstrained Scrsic profile matched to the DOGs using GALFIT (filled grey 
region), DBGs (Distant Blue Galaxies, diagonal hatched blue region) and DRGs (opposite diagonal hatched red region). Right: Distribution 
of i?cff values for DOGs (filled grey regio n), active DRGs (dia gonal hatched blu e region), and quiescent DRGs (opposite diagonal hatched 
red region). DBG and DRG data from ZTrm et al.l l[2003) and lToft etWi | |2007D . 
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TABLE 4 
PSF Subtraction Analysis 





Source Name 


Vnuc 


Inuc 


^^nuc 


gai 


gal 


gai 


ropt 

.^PSF 


j-ir 

■> PSF 


SST24 


J142538.2- 


h351855 


> 28.3 





27±2 


> 26.0 





24.1±0.3 


— 


<0.27 


SST24 


J142622.0- 


h345249 


— 


> 27.9 


27±1 


— 


24.2±0.1 


23.6±0.2 


<0.01 


<0.15 


SST24 


J142626.4- 


h344731 


> 27.6 


— 


26.8±0.3 


> 26.4 


— 


23.7±0.2 


— 


<0.07 


SST24 


J142644.3- 


h333051 


> 28.1 


— 


22.3±0.1 


26.1±0.6 


— 


23.4±0.1 


<0.06 


0.73±0.07 


SST24 


J142645.7- 


h351901 


> 27.7 


— 


27±1 


> 26.2 


— 


22.9±0.1 


— 


<0.08 


SST24 


J142648.9- 


h332927 


— 


> 28.1 


25.2±0.3 


— 


24.9±0.4 


23.3±0.2 


<0.03 


0.15±0.06 


SST24 


J142653.2- 


h330220 


— 


> 28.3 


25.8±0.2 


— 


25.0±0.4 


22.8±0.2 


<0.01 


0.06±0.03 


SST24 


J142804.1- 


h332135 


> 27.7 


— 


26.1±0.5 


> 26.6 


— 


> 25.2 


— 


— 


SST24 


J142924.8- 


h353320 


> 27.9 


— 


25.9±0.5 


> 26.1 


— 


> 25.3 


— 


— 


SST24 


J142958.3- 


h322615 


> 28.2 
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Fig. 6. — Color-magnitude diagram for DOGs, broken down into an extended component (unconstrained Sersic profile) and an unresolved 
nuclear component (from TinyTim PSF). Left: V — H vs. H for DOGs observed by WFPC2 (black squares). The extended component of 
each DOG is shown with a blue circle and the point source with a red star. Detections (Lower limits) are plotted with filled (open) symbols. 
Open green triangles show the median and Icr dispersion in colors of high-z galaxies from the HDF-N (Papovich, personal communication 
and HDF-S l|Labbe et al.|[2003) . The subset of these galaxies qualifying as DRGs are shown with a filled orange triangle. Right: I — H vs. 
H for DOGs observed by ACS. Symbols are same as in left panel. 
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reliable results (|Lotz et al.ll2004f) . In the rest-UV, where 
the DOGs are very faint, none of the 22 WFPC2 im- 
ages and only 6 out of 9 ACS images have the per-pixel 
S/N necessary to compute Rp, G, M20, and C. In the 
rest-optical, however, DOGs are much brighter and 23 
out of 31 NICMOS images have sufficient S/N. Table [5] 
presents the visual and non-parametric measures of DOG 
morphologies, including the per-pixel S/N, Rp, G, M20, 
and C values for the ACS/WFPC2 and NICMOS images. 

In the left panel of Figure [3 we plot G as a function 
of M20 as measured in the rest-UV for DOGs, a field 
galaxy sample, a,nd simulated r^^^ bulges and pure ex- 
ponential disks (|Lotz et al.ll2OO60 . None of the DOGs 
fall within the pure exponential disk- or r^/* bulge- 
dominated regime. The field galaxy population is com- 
posed of sources identified within the ACS FOVs and is 
selected to span the same magnitude range as the DOGs 
in our sample. We use our NDWFS data to apply color 
cuts in Bw — R and R—I space in order to remove objects 
with colors typical oi z < 0.7 sources. The morphologies 
of this field galaxy sample are represented with gray con- 
tours. Four out of six DOGs lie in the lower left corner 
of the plot, with low G and high M20 values indicating 
irregular, diffuse morphologies. In general, low G and 
high M20 values are indicative of dust-enshrouded stellar 
populations, where obscuration by dust causes a galaxy 
to app ear very clumpy wi th flux distributed among many 
pixels ()Lotz et al.ll2008bf) . 

One object (SST24 J143523. 9-1-330706, panel 27 in 
Fig. [3]) has a higher G and lower M20 value than nearly 
all of the field galaxies. This is the same object that 
shows a stronger point-source contribution in the rest- 
UV than the rest-optical. Visual inspection of this ob- 
ject's cutout image reveals an extended feature fading 
towards the southeast that is present in both the rest- 
UV and rest-optical. It appears that the central activity 
in this source is not quite as obscured as in other DOGs, 
but it is not yet clear why this is the case. 

The right panel in Figure [7] shows G and M20 values as 
measured in the rest-optical for DOGs as well as LBGs 
in the HDF-N an d a sample of local {z < 0.1) ULIRGs 
(|Lotz et al.l 120041) . DOGs shift to more typical morpho- 
logical parameters in the rest-optical compared to the 
rest-UV, but they are offset from the parameter space 
occupied by LBGs and local ULIRGs. The median G 
and M20 values for DOGs are 0.49 and —1.24, respec- 
tively, while for LBGs they are 0.63 and —1.6 and for 
local ULIRGs they are 0.59 and —1.5. Part of the differ- 
ence in G and M20 compared to LBGs may be that LBGs 
are more compact, and hence less resolved. The offset 
to lower G values in DOGs compared to local ULIRGs 
is remarkable and indicates that either different mecha- 
nisms are involved in creating these two populations, or 
they represent different stages in the evolution of massive 
galaxies. We note that the ULIRG sample has compa- 
rable S/N as the DOGs studied here, and that the rest- 
frame wavelength of both samples is similar (~7000 A vs. 
~5300 A, respectively). While there is a greater rela- 
tive difference in the spatial resolution of the two sam- 
ples (^0.2 kpc pi x~^ vs. (saO.6 kpc pix~^, respectively), 
iLotz et all ()2004[ ) found that systematic offsets at these 
resolutions should be on the order of 20% or less for G 
and M20 and less than 10% for G. Therefore, the dif- 



ference in G cannot be explained by spatial resolution 
effects alone. 

As a qualitative consistency check, we examined 
i?-band imag e s of a sample of 56 ULIRGs from 
iMurphv et all (1996) and determined that 20 (35%) 
have double nuclei with nuclear separations larger than 
2.3 kpc, approximately the resolution limit of our NIC2 
images. In comparison, only one DOG has two well- 
detected, distinct nuclei and three or four have low S/N 
components separated by 0'.'5 (Ri4 kpc), implying that 
at most 16% of the DOGs in our sample have multiple 
nuclei with separations larger than 2.3 kpc. This result 
is qualitatively consistent with the differences seen in G 
between local ULIRGs and DOGs. An important caveat 
with this analysis is that our sample of DOGs is domi- 
nated by power-law sources, while the ULIRG sample has 
a variety of rest-frame NIR SED shapes. For reference, 
we measured the G and A/20 values of a well-detected, 
non-DOG point source (S/N per pixel of 18) in one of 
our NIC2 images, and found values of 0.62 and —1.7, re- 
spectively. The DOG whose morphology is dominated 
by a point source (see panel 31 in Fig. [3]) has a lower G 
value (0.56) but almost the same M20 (—1.8). This may 
be an indication that this DOG contains an underlying 
extended component, but the data are not conclusive. 

5. DISCUSSION 
5.1. Dust and Stellar Mass Estimates 

Here we estimate some of the intrinsic properties of 
DOGs, including lower limits on their reddening (Ay), 
dust and gas masses, and stellar masses. To do this, 
we use Si mple Stellar Popula t ion (S SP) template SEDs 
from the iBruzual fc Chariot! (|2003l ) population synthe- 
sis library with ages spaced logarithmically from 10 Myr 
up to 1 Gy r, as w ell as the median QSO template from 
lElvis et all (|1994f ). All models used here have solar 
metallic ity, a C habrier IMF over the mass range 0.1- 
IOOM0 (|Chabri er 2003) , and use th e Padova 1994 evolu- 
tionary tracks (jCirardi et al.lll996[). The redd ening law 
used is a combination of that from lCalzetti et al. (200(|) 
and longer wavelength estimates from Draind ()2003) , and 
assumes the case of a dust screen in front of the emitting 
source in order to derive a &m lower limit on Ay- 

For each DOG, we estimate Ay needed as a function of 
age by determining the amount of extinction necessary 
to redden the given SSP template such that it repro- 
duces the observed V — H 01 I — H color. The process 
is illustrated in Figure [51 Each panel shows the colors 
of DOGs in our sample as a function of redshift. Blue 
circles represent the extended component and red stars 
show the point source component of each DOG, as de- 
scribed in section 14.2.11 Dotted lines show the expected 
colors of the SSP templates for varying amounts of ex- 
tinction. Even with no extinction {Ay = 0), the oldest 
SSP templates are too red to reproduce the colors exhib- 
ited by DOGs. The QSO templates, on the other hand, 
require large Ay values in order to match the DOG nu- 
clear colors. 

We use the relation from iBohlin et ahl (|1978f ) to con- 
vert Ay to the total column density of hydrogen atoms 
and molecules, A^h- For the 100 Myr SSP, the column 
densities range from 3 x 10^" — 6 x 10^^ cm^^, with a 
median A'h of 3 x 10^^ cm^^. If we assume the dust is 
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TABLE 5 

Non-parametric Morphological Classifications 
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^ Mode of visual classification. Italics indicate multiple users disagreed with the mode. Negative C value indicates r20 was too small 
to be measured accurately. 




Fig. 7. — Gini coefficient vs. M20- Left: Morphological measures from ACS/F814W images of DOGs (blue squares) and field galaxies 
(grey solid contours, see text). Representative error bars in the lower right corner include uncertainties due both to low S/N and low spatial 
resolution and are estimated using a method similar to that in Lotz ct al. (2004). Filled symbols have greater than 20% point-source 
contribution. Top and bottom dotted boxes show where simulated face-on bulges and disks lie, respectively IILotz et al.ll200(3) . Right: Same 
plot but showing r esults from NIC2 /F160W images of DOGs (red circles), HDF-N LBGs (filled green star) and a sample of local ULIRGs 
(black diamonds) l ILotz et al.|[200g) . 
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Fig. 8. — V — H (top row) and I — H (bottom row) as a function of spectroscopic redsliift for eacli DOG. Tiie first three columns 
show the colors of the extended component, while the fourth column shows the colors of the unresolved component (filled symbols are 
detections, open symbols are lower limits). Filled orange triangle s represent DRGs in the H DF-S. Dotted lines trace the evolution of colors 
with redshift of reddened simple stellar population models from lBruzual fc Chariot 1 12003 ) with solar met allicty, a Chabrie r IMF and at 
ages of 10 Myr, 100 Myr, and 1 Gyr (three columns on left), as well as of the median QSO template from .Elvis et al] I) 19941 '). 
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distributed in a spherical shell around the source with 
radius equal to the effective radius, then we can place a 
lower limit on the dust mass: 



1 9 

Mdust > -^Mp^H X 47ri?^ff. 

/gd 



(1) 



Here, /gd is the gas-to-dust mass ratio and fip is the mean 
molecular weight of the gas, which we take to be l.Bmp, 
where rup is the mass of a proton. We adopt the average 
gas-to-dust mass ratio of 120 measured for th e nuclear 
regions of local ULIRGs bv lWilson et"aI1 (|200l . We find 
dust mass lower limits ranging from 2 x 10^ — 6 x lO^M©, 
with a median of 9 x 10^ Mq for a 100 Myr SSP. 

A complementary method of estimating the dust mass 
is based on measurements of optically thin sub-mm emis- 
sion. Because sub-mm photometry for the DOGs in our 
sample is currently unavailable, we extrapolate from the 
24/im flux density measurement. We used the Mrk231 
template to determine the extrapolated 850/im flux den- 
sity, ^850 . Using a template SED of a galaxy with colder 
dust such as Ar p220 would inc r ease t he inferred 850/im 
flux. We follow iHughes et all (|1997[ ) and estimate the 
dust mass using: 



(2) 



where is the luminosity distance, Kd is the rest- 
frequency mass absorption coefficient, and T) is the 
value of the modified blackbody function (/3 = 1.5) at 
the rest frequency v and a temper ature T. The appro- 
priate Kd value is interpolated from iDraind (|2003f ). with 
typical values being 5 cm^ g"-'^. There is at least a fac- 
tor of 2 uncertainty in this quantity. We have assumed 
relatively hot dust (Td = 75 K), since we expect AGN 
heating to play an important role in this sample of DOGs 
(a dust temperature of 50 K would increase the inferred 
dust mass by an additional factor of «1.5). Using this 
method, we find dust masses of 8 x 10 - 6 x 10^ Mq, 
with the median dust mass being 1.6 x lO'^ Mq. This is 
a factor of nearly 20 larger than the median dust mass 
inferred from the measurements of Ay. This might be 
expected, given that many of the dust masses based on 
Ay are lower limits, while the dust masses based on the 
24/im emission may be overestimates if Td > 75 K. On 
the other hand, this difference may be suggesting that 
the dust causing the average UV extinction of the ex- 
tended galaxy component is not the same dust that is 
causing the thermal emission. 

We use the SSP templates to estimate the stellar mass 
in each DOG. This is computed by reddening each SSP 
template to match the observed color of the DOG at 
the appropriate redshift. We then scale the redshifted, 
reddened template to match the obse r ved H -band pho- 
tometry. Since the lBruzual fc CharlotI (|2003[ ) models are 
normalized to a stellar mass of 1 Mq, this scaling factor 
represents the stellar mass of the DOG. 

In FigureEl we show the stellar mass of each DOG as a 
function of age as well as the distribution of stellar masses 
assuming a 100 Myr SSP model. As stellar populations 
age, their colors naturally redden, and so less extinction 
is needed to reproduce the observed colors of the DOGs. 
For most DOGs, ages greater than ~300 Myr require Ay 
values less than zero and are unphysical. Meanwhile, at 



younger ages, lower mass-to-light ratios are balanced by 
the need for greater extinction to match the observed col- 
ors. As a result, the inferred stellar masses are relatively 
constant to within a factor of a few for ages less than 
300 Myr. We note that these mass estimates are lower 
limits because (1) the amount of extinction is a lower 
limit, especially when there is no detection in the V- or 
/-band image and (2) our extinction estimate does not 
take into account grey extinction. For an age of 100 Myr, 
the stellar masses range from 2 x 10® — 1 x 10^^ Mq, with 
the median mass being 3.3 x 10^° Mq. 

We use our dust mass estimates inferred from the 
24/im flux den sities and a gas-to-dust mass ratio of 120 
(jWilson et al.| [2008) to obtain an upper limit on the gas 
masses. This leads to gas masses of l-7xlO^''M0, with a 
median gas mass of 2 x 10^". If we assume a closed-box 
model and an exponential star-formation history, then we 
can write Mgas -|- Mdust = Mtotexp(-t/to), where Mtot 
is the sum of the gas, dust, and stellar masses. The stel- 
lar mass is then given by Mgtar — Aftot(l — exp(— i/io)), 
which implies 



tn 



Mt, 



M, 



dust 



(3) 



This quantity represents the fractional lifetime (in 
units of the scale time for our exponential star-formation 
rate assumption) of each DOG. Larger values of t/to in- 
dicate more evolved systems, as more of the gas has been 
converted to stars. Our values of i/^o represent lower lim- 
its on the actual values because our stellar masses are un- 
derestimated and our gas masses are based on our 24/im 
flux densities, which likely overestimates the true mass 
of gas. We find t/to lower limits ranging from 0.02 to 
3.3, with a median lower limit of 0.9. This result implies 
that in half of the DOGs in our sample, at least 90% of 
one exponential timescale's worth of star-formation has 
occurred. The "oldest" DOG in our sample has gone 
through more than three exponential timescales of evo- 
lution. However, we caution that the "youngest" DOGs 
from this line of analysis are uniformly associated with 
sources where the dust extinction is most likely under- 
estimated, thereby causing an additional underestimate 
in the stellar mass and the associated t/to value. In Ta- 
ble [6l we present the dust and stellar masses derived in 
this section, as well as our measure of the lower limit on 
the fractional lifetime, t/tg. 

5.2. Comparison to other high redshift galaxy 
populations 

It is important to understand how DOGs are related to 
other populations of high-z galaxies that have been stud- 
ied in the literature. Here we compare the morphological 
properties of the DOGs with some of these high-redshift 
galaxy populations and find that DOG morphologies are 
distinct from the bulk of LBGs and quiescent high-z 
galaxies, but are similar to SMGs as well as active DRGs 
and the extreme subset of faint, diffuse LBGs. 

5.2.1. Sub-mm Galaxies 

SMGs are a particularly interesting population of 
galaxies to compare with DOGs. First identified by blind 
sub-mm surveys with the Submillimetre Comm on User 
Bolometer Array fSCUBA iHoUand et al.|[l999l ). SMGs 
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Fig. 9. — Left: Stellar mass as a function of SSP age for the DOGs. Right: Distribution of stellar masses at an age of 100 Myr. 
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may represent an important, short-lived, and very ac- 
tive phase in the evolution of the most massive galaxies 
Their redshift distribution, number density, an d clus- 
terin g properties are s imilar to DOGs ( [Chapma n et al 



20051 iDev et all 120081 iBlain et alJ l2004iriBrodwin et al 



2008| ). However, a sample of DOGs detected at TO/um 



or 160/zm by Spitzer/MlPS tend to show warmer colors 
(i.e., smaller 70/24/im or 160/24/im flux density ratios) 
compared to SMGs (Tyler et al., submitted). One spec- 
ulative scenario that may serve as a possible explanation 
for this behavior is that these two galaxy populations are 
linked in an evolutionary sense: SMGs represent a cold 
dust, star- formation dominated stage in the formation of 
massive galaxies that may precede the DOG phase, when 
the feedback from the growth of a central black hole has 
heated the surrounding gas and dust, thereby quench- 
ing star formation and shifting the peak of the SED to 
shorter wavelengths. 

If this scenario is correct, then we expect to see 
major mergers dominate the morphologies of SMGs, 
while DOGs should show more relaxed morphologies 
typical of the f inal m erger stage before the remnant. 
IConselice et al.l ()2003l ) analyzed Space Telescope Imag- 
ing Spectrograph (STIS) rest-frame UV data of a sample 
of 11 SMGs at z ~ 2 — 3 using the CAS system, and found 
evidence suggesting a major merger fraction of between 
40% and 80%. Although we do not have the S/N in our 
images to measure A reliably (and thereby determine a 
major merger fraction in a similar manner), the low G 
and high M20 values we have measured imply diffuse, ir- 
regular systems where the light is spread into multiple 
components rather than 2 separate components. If the 
DOGs were predominantly major mergers, we would ex- 
pect our sample to have higher G values. Instead, the 
low G values we find suggest that we may be looking 
at the dusty remnant of a major merger, where there 
are many highly obscured components near each other. 
However, we caution that dust can have a strong effect 
on the measured G value in the rest-UV, such that even 
major mergers might yield lower G values. 

The rest-UV morphologies of SMGs have also been an- 
alyzed in the GOODS-N field, where a SCUBA super- 
map exists and has been used to identify robust sub-mm 
detections (jBorys et al.ll2003[ ). A sample of 12 sources 
in the redshift r ange 1.7-4.0 (comp arable to the DOGs) 
were studied bv lPope et al.l (j2005[ ). who computed con- 
centration and asymmetry values, finding C to be in the 
range 2-3.3 and A to be dominated by noise, with the 
exception of two objects (one is very compact and the 
other is clearly asymmetric). The comments associated 
with many of these sources are "faint" and "diffuse" , sug- 
gesting that there is large-scale dust obscuration in these 
systems. This is qualitatively similar to what is seen in 
many of the DOGs, suggesting that there is some overlap 
between the two samples. 

SMGs have not yet been characterized in terms of G or 
Af20i so direct comparisons based on these quantities are 
not possible at this time. However, we can compare the 
sizes of these systems directly via the Petrosian radius. 
In the rest-UV, the DOGs range in size from Rp ~ 0.5" 
to 1.5", while SMGs range in size from 0.5 - 2.5". Indeed, 
a two-sided KS test reveals that there is only a 5% chance 
that they are drawn from the same parent distribution. 
This suggests that, while there are similarities between 



SMGs and DOGs, SMGs tend to be slightly larger than 
DOGs. This is consistent with the major merger hypoth- 
esis in which DOGs are in a more evolved state where 
dynamical friction has caused individual components to 
fall towards the center of mass. However, we caution 
that this result in itself does not provide evidence for 
DOGs originating from major mergers. Objects mov- 
ing at ~100 km will traverse 8 kpc in <100 Myr. 
Simulations of major mergers predict a phase of intense 
star formation and central black hole growth that lasts of 
order this timescale, indicating that the size differences 
are at least consisten t with the scenario outlined above 
(jHopkins et al.ll2007f ). 

Finally, the ratio of the stellar mass to the gas mass 
holds potential for comparing the evolutionary states of 
SMGs and DOGs. The ideal comparison study would 
include statistically significant samples of both popula- 
tions of galaxies for a range of observed properties such as 
24/im emission, bolometric luminosity, space density, etc. 
Unfortunately, such samples do not currently exist — 
mainly due to a combination of the limitations of current 
instrumentation and the fact that these are relatively 
recently discovered populations of galaxies. While CO 
linewidths and emission strengths have found gas mass 
estimates for a handful of SMGs, no such measurements 
have been published for DOGs. The gas mass estimates 
for SMGs are tvpicallv ^ 5x10^° Mq ICreve et al.l [20051 : 
iTacconi et al[|2006l |2008[) . A number of efforts have been 
directed at determining the stellar masses of SMGs using 
SED-fitting algorithms. Average stellar mass val ues are 
in the range 3-6x10" M© (Borys et al. 2005; Dv e et al.l 
l2008f ). However, when we employ our method of deter- 
mining the stellar mass (in this case using the R—K color 
to determine the optimal Ay value for a given SSP and 
age) using the photometry presented in those papers, we 
find average stellar masses of ^7x10^*^ Mq. These esti- 
mates increase by a factor of ^2 if instead of a Chabrier 
IMF we use a Salpeter IMF (as was done by the previ- 
ous authors for SMGs). However, this still leaves us a 
factor of k2 short of the mass estimates provided in the 
papers described above. In order to compare DOG stel- 
lar masses with SMGs consistently, we adopt the lower 
stellar mass values that we derive for SMGs. In this case, 
the median t/to values for SMGs becomes ^lA, which is 
much closer to the median lower limit value of 0.9 found 
for the DOGs in section [5TT1 

The large uncertainty inherent in the process of es- 
timating dust and gas masses based on 24/Ltm photom- 
etry or rest-frame optical Ay measurements currently 
prevents a strong conclusion being made regarding the 
evolutionary status using this line of analysis. However, 
the morphological evidence is suggestive — although not 
conclusive — of an evolutionary link between the two 
populations with SMGs serving as the less evolved pre- 
cursor to the DOG phase. 

5.2.2. Star-forming Galaxies 

A number of selection criteria have been used to 
identify normal star-forming galaxi es at high redshift . 
Two of these are t he LBG dro pout (iMadau et al.lll996l : 
ISteidel et all I1996D and BzK f Daddi et al.l l2004f ) tech- 
niques. A direct comparison to our work can be made 
with a sample of LBGs an d emission hne ga laxies in the 
GOODS-N field studied bv lLotz et all ((200I) . These au- 
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thors compared G, M20, and C values between their sam- 
ple of 82 z ~ 4 LBGs and 55 z 1.5 emission line galax- 
ies. In the LBG sample, they found a major-merger frac- 
tion of ~ 10 - 25% (defined by M20 > -1.1) and a bulge- 
dominated fraction of - 30% (G > 0.55, M20 < -1.6). 
The remainder of the LBGs had G and M20 values larger 
than what is typical for normal galaxies, suggesting ac- 
tive star-formation or a recent merger event. The low- 2: 
emission-line sample showed a similar major merger frac- 
tion but fewer bulge-dominated systems. It is remark- 
able then, that so few of the DOGs have G and M20 
values typical of bulge-dominated systems, even in the 
rest-optical, despite their luminosity. Furthermore, four 
out of six DOGs with measureable morphologies in the 
rest-UV have high M20 and low G values that are typical 
of dusty, irregular systems. This may be an indication 
of kpc scale dust obscuration, which can bias the G and 
M20 values away from the bulge-dominated regime. 

A morphological study of LBGs in GOODS-N by 
iRavindranath et al.l (|200 6^ found axial ratios skewed to- 
wards lower values for galaxies at 2; > 3, suggesting high- 
z LBGs are dominated by edge-on morphologies. In con- 
trast, only one DOG has an axial ratio less than 0.35, 
indicating that if these sources are disk galaxies, then 
some selection mechanism must be in place that favors 
observing DOGs in face-on orientations. Meanwhile, re- 
sults from numerical simulations of galaxy mergers indi- 
cate that remnants end up with axial ratios between 0.5 
and 1.0, depending on the viewing angle (jNovak et al.l 
12006( 1 ■ 21 DOGs satisfy this axial ratio criterion, but the 
median value in our dataset is «0.5. This suggests that 
either DOGs represent a phase prior to the final remnant 
stage or th ey are form ed by some other process. 

Recently, iLaw et al.l ([2007) have used GOODS data to 
analyze morphologies of 216 LBGs and compare them 
with other high-z galaxy populations. They found sig- 
nificant overlap between the LBGs and BzKs, indicating 
that the optical and NIR selection criteria are identify- 
ing similar galaxies. While these authors performed a 
non-parametric morphological analysis, direct compari- 
son between our work and theirs is difficult because (a) 
they assign pixels to each galaxy based on an isophotal 
surface brightness criterion rather than the elliptical Pet- 
rosian radius as we have done and (b) they create their 
own parameter to describe the multiplicity (^E") of each 
galaxy, rather than using M20. Nevertheless, there are 
some ap parent diffe r ences between the DOGs and LBGs 
from the ILaw et all (|2007t l study. While the LBGs span 
the full range of G values in the rest-UV, the DOGs tend 
to be low G objects. Furthermore, though LBGs span 
a wide range in G, they are preferentially found to have 
low ^E* values, implying a small number of distinct compo- 
nents. On the other hand, the DOGs have high A/20 val- 
ues, s uggesting multi-co mponent structure is common- 
place. ILaw et al.l ()2007( ) note a correlation in their plot 
of G as a function of 4", in the sense that objects with 
many components (large tend to be fainter and more 
nebulous (low G). DOGs resemble this extreme subset 
of faint diffuse LBGs, but appear highly morphologically 
distinct from the vast majority of the LBG population. 

5.2.3. Passively Evolving Galaxies 

As mentioned above, the BzK method can be used 
to identify high- 2: passively evolving galaxies. This pho- 



tometric color cut has been used in the Hubble Ultra 
Deep Field (UDF) to generate a samp le of seven luminou s 
early-type galaxies at 2; = 1.39-2.47 ()Daddi et al.ll2005D . 
These authors studied the i and z band morphologies 
of all seven objects with both parametric (Sersic profile 
fitting) and non-parameteric (concentration and asym- 
metry) methods. They found fairly large Sersic indices 
(n ~ 3) and small effective radii (reff < 1 kpc), typical 
of E/SO galaxies. In contrast, the best-fit Sersic profile 
for DOGs has smaller n values more typical of expo- 
nential disks (median n = 0.9) and larger effective radii 
(i?cfr '^1 — 6 kpc). Moreover, the passive BzK galaxies 
have G > 2.6 and A < 0.2, consistent with early-type 
systems. In the rest-optical, DOGs tend to show lower 
G values (S/N is not sufficient to measure A), consis- 
tent with an exponential profile. Along with the low G 
and high M20 values that are measured for the DOGs, 
these morphology results suggest that DOGs and pas- 
sively evolving high-z galaxies are distinct populations, 
either because they represent different stages of evolution 
or because they have different formation mechanisms. 

5.2.4. Distant Red Galaxies 

Another population of high-z galaxies is the so-called 
Distant Red Galaxies (DRGs). Identified via deep NIR 
imaging, these objects were fi rst postulated to be the 
reddened descendents of LBGs (jFranx et al.ll2003l ). Sub- 
sequent studies of DRGs in the Extended Groth Strip 
(EGS) show a wide variety of shapes, with 57% ap- 
pearing visually as elliptical/compact, 7% as edge-on 
disks, and the remai nder as peculiar/irregular galaxies 
(jConsehce et al.ll2007l) . The low redshift DRGs (z < 1.4) 
have GAS values typical of nearby normal galaxies. The 
higher z DRGs visually classified as elliptical/compact 
have higher G values, similar to what is seen locally in 
massive elliptica ls and in the BzK samples. Meanwhile, 
ILaw et all (|2007[) examined DRGs in the GOODS-N field 
that did not overlap with the BX or BM LBG color cri- 
teria and found that this population of galaxies was sub- 
stantially fainter and more diffuse than either the star- 
forming BzKs or the LBGs. They note that this is the 
behavior one expects from dusty, IR-bright galaxies. The 
faint, diffuse nature of these objects is reminiscent of the 
DOGs, and it is possible that there is significant overlap 
between these two populations. 

Previous work using rest- frame UV-NIR SEDs has sep- 
arated actively star- forming DRGs (sDRGs) from quies- 
cent ones (qDRGs) (jZirm et al.ll2007t iToft et al.ll20071 
Examination of the morphological differences between 
these two populations has revealed a correlation between 
size and star formation activity, in the sense that qDRGs 
are all very small (i?cff ^ 1 kpc), while sDRGs span a 
larger range in size (i?cff ~ 1 — 10 kpc). As is shown in 
the right panel of Figure [5l DOGs appear very similar 
to sDRGs in terms of their sizes. This is consistent with 
the qualitative similarity between sDRGs and DOGs de- 
scribed in the preceding paragraph and suggests there is 
extensive overlap between these two populations. 

5.3. Implications for the Evolution of the Most 
Massive Galaxies 

In the local universe, there has been evidence for some 
time that warm dust-dominated ULIRGs may represent 
a transition stage between cold ULIRGs and optically 
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luminous quasars (jSanders et al.lll988bD . If this scenario 
holds at high redshift, then there is a natural explana- 
tion for the observations based on the selection criteria 
alone: objects selected at long wavelengths (i.e., SMGs) 
are preferentially cold-dust dominated systems and rep- 
resent the 'cold ULIRG' phase, whereas objects selected 
at 24/im (i.e., DOGs) are dominated by warmer dust 
and represent the transition phase en route to the op- 
tically luminous quasar. As time progresses and the 
quasar fades in luminosity, the compact, quiescent, el- 
liptical galaxy remnant becomes visible (i.e., quiescent 
BzKs and DRGs). 

Again referring to the local universe for guidance, if the 
triggering mechanism for this activity is a major merger 
(jSanders et al.l [l988a[ ). then we should expect to see a 
trend in relaxation and size, where the initial stage shows 
the largest sizes and least relaxation and the end prod- 
uct is a relaxed, compact system. This picture is ap- 
parently consistent with our data, as DOGs tend to be 
smaller than SMGs, but larger than quiescent DRGs or 
BzK galaxies. Furthermore, SMGs frequently exhibit 
signs of major merger activity, whereas passively evolv- 
ing systems at high-z are very compact with large Sersic 
indices. DOGs appear to be intermediate stage objects 
that typically do not show signs of major mergers, but 
nonetheless have morphologies indicating they are more 
dynamically relaxed than SMGs but less so than the qui- 
escent systems. 

It is important to emphasize that while our morpholog- 
ical results are consistent with the hypothesis that DOGs 
act as a transition phase in the process of creating a mas- 
sive galaxy via a major merger, the morphological infor- 
mation currently available is not sufficient to exclude the 
possibility that DOGs are created by some other process 
such as minor merging (for example, minor mergers have 
the potential to increase size temporarily), or are simply 
dusty galaxies hosting a powerful, obscured AGN. 

Our analysis of the stellar, dust, and gas masses of 
DOGs currently does not provide compelling evidence to 
place them within an evolutionary scheme with respect 
to other massive proto-galaxy candidates such as SMGs. 
Additional data are needed before conclusive statements 
can be made based on mass estimates such as these. 

6. CONCLUSIONS 

We have analyzed the morphologies of 31 Dust Ob- 
scured Galaxies (DOGs) at z « 2 from the Bootes field 
using data from TJST ACS/WFPC2 and NICMOS. Our 
findings are summarized below. 

1. Although these sources were selected to have mid- 
IR signatures of AGN, we detect spatially resolved 
emission at rest-frame UV and/or rest- frame opti- 
cal wavelengths for all but one of the 31 targets. 

2. Using a three component model in GALFIT (sky 
+ PSF -|- Sersic profile), we measure significant un- 
resolved components in 28 out of 31 DOGs in the 
rest-optical, and the median point-source fraction 
is 0.13. Only 10 DOGs have measureable unre- 
solved components in the rest-UV. 

3. The median Sersic index is 0.9, indicating that 
disk-like profiles are preferred to bulge-like ones. 



On the other hand, very few DOG extended com- 
ponents have small axial ratios, indicating that if 
DOGs are predominantly a population of normal, 
disk-like galaxies (with an obscured AGN produc- 
ing the 24/im flux), then some selection mecha- 
nism(s) must be in place that favors face-on rather 
than edge-on orientations. 

4. DOGs in our sample have effective radii of 1-5 kpc, 
which places them between SMGs and quiescent 
DRGs or BzK galaxies. If DOGs are formed by a 
major merger, this trend in sizes is consistent with 
them acting as a transition stage in the evolution 
of massive galaxies. If DOG activity is triggered 
by some other process, such as a minor merger or a 
dusty AGN in a normal galaxy, then interpretation 
of this size trend is not as clear. 

5. In the rest-optical, DOGs have lower G values than 
local ULIRGs (median values of 0.49 and 0.59, re- 
spectively). This might be expected if DOGs repre- 
sent a subsequent stage in the merging process (just 
before coalescence), but might also be expected if 
the galaxies are not disturbed by a major merger. 

6. Simple stellar population modeling reveals that old 
(>300 Myr) single-burst stellar populations are 
redder than most DOGs and thus ruled out. If 
100 Myr old SSPs are appropriate, then DOGs re- 
quire substantial amounts of extinction to produce 
the observed red colors, with Ay — 0.2 — 3. This 
provides a lower bound on the median dust mass of 
10'' Mq. An upper bound is obtained by extrapo- 
lating the 24/im fiux density to 850/im and is found 
to have a median value of 1.5 x 10^ Mq. We find 
a median stellar mass lower limit of 3 x 10^*^ Mq 
which is relatively insensitive to age to within a 
factor of a few. 
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